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^ , Abstract. We test the canonical BCS wave functions for fixed number of electrons for the 

Oh' attractive Hubbard model. We present results in one dimension for various chain lengths, 

^ , electron densities, and coupling strengths. The ground-state energy and energy gap to 

the first excited state are compared with the exact solutions obtained by the Bethe ansatz 
C^ , as well as with the results from the conventional grand-canonical BCS approximations. 

While the canonical and grand-canonical BCS results are both in very good agreement 
with the exact results for the ground state energies, improvements due to conserving the 
electron number in finite systems are manifest in the energy gap. As the system size 
is increased, the canonical results converge to the grand-canonical ones. The "parity" 
^ ' effect that arises from the number parity (even or odd) of electrons are studied with our 

canonical scheme. 



o 



> 

•^ \ In recent experiments on nanoscale superconducting Al particles (e.g., Ref. [1]), 

^ ' the number of electrons in a particle was fixed by charging effects and the excitation 

CS ■ spectrum has been measured by tunneling of a single electron. They have observed the 

number "parity" effect, the difference in the excitation spectra that arises due to the 
Q\ \ number of electrons in a particle being even or odd. For these systems, the validity of 

the BCS theory that utilizes the grand canonical ensemble is questionable. We address 
this issue by using the attractive Hubbard model, which serves as the minimal model 
to describe s-wave superconductivity. We test the canonical BCS wave function that 
'^ \ has a fixed number of electrons for this system. While our formulation can be applied 

for larger systems up to three dimensions, we present results for one dimension, for 
O ■ various system sizes, electron densities and coupling strengths. Exact solutions are 

available in one dimension via Bethe Ansatz techniques [2], and the grand canonical 
\rl \ BCS solutions have recently been evaluated for large system sizes [3]. 

;h ' The attractive Hubbard Hamiltonian is given by 



c^ 



c^ 



H = -J2ts (a^^Oi,, + h.c.) - \U\Y,ni^nii (la) 

^5 i 
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= II ek aL^ka -^Y. «kT«-k+u«-k'+uak'T , (lb) 

kcr ^^ kk'l 
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where a]^ (oio-) creates (annihilates) an electron with spin a at site i and riia is the 
number operator, Uia = di^ciia {i is the index for the primitive vector Rj). The t^ is the 
hopping rate of electrons from one site to a neighbouring site R^ away. In this work, 
we include the nearest neighbours only, and t = ts- The \U\ is the coupling strength 
between electrons on the same site, and we have explicitly included the fact that the 
interaction is attractive. In Eq.(lb), we have Fourier transformed the Hamiltonian 
with periodic boundary conditions for A^ sites in each dimension. The kinetic energy 
is given by 

ek = -2^t5Cosk-R5, (2) 

5 

where R^ is the coordinate vector connecting sites i to i + 6. 

We perform variational calculations using the component of the BCS wave function 
that has a given number of pairs z/ and thus conserves the number of electrons Ng = 
2v. The BCS wave function is a superposition of pair states with all the possible 
numbers of pairs {z/}: 

oo 

|BCS)gc = n(«k + t;k4^alkj |0) = E 1^-) ' (3) 

k v=0 

where |0) denotes the vacuum state, and |\E'o) = |0). The z/-pair component \^u) can 
be obtained by rearranging |BCS)gc into a power series of the pair creation operator 
^kt'^-ki ^^*i c^^ be written as 

l^-) = ^n(E^k,aLT«Ui)|0)> (4) 

'^- i=i k, 

where kj ^ k^ for all i ^ j, and we have defined g]^- = (llk'^k) '^ki/'^kr An 
alternative way of formulating the problem is to express the above wave function as a 
particle- number projection of |BCS)gc [4^6]. The wave function with an odd number 
of electrons Nf, = 21/ + 1 is defined by 

i=l ki 

with kj 7^ kj for all i ^ j- We construct the energy gap by the definition 

Aiv. = 1{En.~i - 2Em^ + E^,+i) . (6) 

In Fig. 1 we show the ground state energy as a function of the electron density 
n = Ne/N for A^ = 16, for \U\/t = 10 and 4. The improvements of the canonical over 
the grand canonical energies can be seen, while the former with even N^s converge to 
the latter for smaller n and larger \U\, and both results are in very good agreement 
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FIGURE 1. Ground state energy as a function of the electron density for N — 16. 
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FIGURE 2. Energy gap as a function of the electron density for TV = 32 and \U\/t — 2. 



with the exact ones for small n. The difference between the energies with even and 
odd A^e's is apparent for large \U\. 

The energy gap (scaled by the kinetic energy parameter t) is shown in Fig. 2 
as a function of the electron density, for A^ = 32 and \U\/t = 2. The canonical 
method improves the grand canonical one significantly for small system size and 
weak coupling. As the system size or the coupling strength is increased, the canonical 
results converge towards the grand canonical ones. For \U\/t = 2, N = 32 is large 
enough that the canonical gaps are closer to the grand canonical ones than to the exact 
solutions, and the canonical and grand canonical gaps are more or less converged for 
small n. The gap for odd N^ is negative. Its magnitude is different from the even A'e 
gap for larger density and weak coupling. Also in this limit, the even gap oscillates as a 
function of the density. Both of these effects result from the quantized (unperturbed) 
energy levels due to finite size, and become negligible for strong coupling and for large 
system size (small density). 
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FIGURE 3. Occupation probability for different k states as a function of the coupling strength for 
N = 64a.ndiy = 32. 

Finally, in Fig. 3, we show the occupation probability of different k states as a 
function of the coupling strength. In the strong coupling limit, all the unperturbed 
states are more or less equally occupied (in Fig. 3, this limit is not reached yet). 
As the coupling strength goes to zero, the distribution function for free electrons is 
recovered. Though the figure shown has been produced by the canonical scheme, the 
grand canonical approximation also yields the same probability distribution for such 
a large system size. 

In conclusion, the canonical method yields improvements in the energy gap for 
small system size and weak coupling. Further study is being made for comparison 
with the grand canonical scheme that conserves the number parity [7]. 
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